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Electronic devices lose efficacy due to quantum effect when the line-width of gate decreases to sub- 10 nm. 
Spintronics overcome this bottleneck and logic gates are building blocks of integrated circuits. Thus, it is 
essential to control electronic transport of opposite spins for designing a spintronic logic gate, and 
spin-selective semiconductors are natural candidates such as zigzag graphene nanoribbons (ZGNR) whose 
edges are ferromagnetically ordered and antiferromagnetically coupled with each other. Moreover, it is 
necessary to sandwich ZGNR between two ferromagnetic electrodes for making a spintronic logic gate and 
also necessary to apply magnetic field to change the spin orientation for modulating the spin transport. By 
first principle calculations, we propose a method to manipulate the spin transport in graphene nanoribbons 
with electric field only, instead of magnetic field. We find that metal gates with specific bias nearby edges of 
ZGNR build up an in-plane inhomogeneous electric field which modulates the spin transport by localizing 
the spin density in device. The specific manipulation of spin transport we have proposed doesn't need 
spin-charge conversion for output and suggests a possible base for designing spintronic integrated circuit in 
atomic scale. 

More than thirty years ago, it has been proposed that quantum effect can be used for calculations' and lots 
of experimental and theoretical studies on spintronics emerged". For spintronic logic gates, ZGNR can 
be used as both devices and leads due to its spin-selectivity"*"^ and tunability by electric field*" However, 
it limits the feasibility that ferromagnetic electrodes and magnetic field are used in spintronic logic gates to 
manipulate the spin transport*' Thus, it remains unknown to manipulate the spin transport for logic operations 
('AND', 'OR', 'NOT') with applied electric field directly. Here we show a method of utilizing only electric field to 
tune spin transport in ZGNR, and design a prototype of spintronic logic gate by first-principle calculations. 

In experiments, graphene nanoflakes (GNF) have been fabricated and observed'"' For GNF samples, zigzag- 
shaped edges are more favorable than armchair-shaped edges in energy due to less dangling bonds. Selected 
graphene flakes are further processed into multi- terminal devices"'". An ideal zigzag-shaped edge of ribbon-like 
GNF is ferromagnetically ordered". In real GNF, unordered defects at edges localize the magnetism of edge states, 
although there are lots of Dirac fermions in bulk^" Both good quality with less defects and accurate measure- 
ments are needed to seek the ferromagnetism at edges in GNF. It is reported recently that spin-polarized transport 
at edges in ribbon-like GNF have been observed^''. The experiment (Ref. 23) has shown a new kind of one- 
dimensional electronic system with a tunable band gap and a spatially separated spin texture due to breaking the 
symmetry of planar spin-rotations under a very large magnetic field angled with respect to the graphene plane. In 
principle, the spin texture of ZGNR can be modulated not only by the symmetry breaking in magnetic field but 
also by a lot of other factors such as external electric field"* magnetic impurities"''' ''\ substrate adsorption'", and 
other structure defects^^ ""*. Here we propose a prototype of quantum logic gates based on the following mech- 
anism: Spin polarized states in the vicinity of Fermi level (£p) are localized at edges of GNF, so an external electric 
field localized nearby the edge can be used to tune the distribution of spin density and tune the spin transport 
correspondingly. 

A graphene nanoribbon with four zigzag chains (4ZGNR) is used to construct the logic gate whose size is much 
smaller than the gate line-width of current CPUs. At each edge of 4ZGNR, two metal gates are established and a 
dielectric region is inserted to divide the external potential at the edge into three parts (Fig. la). For our model, the 
dielectric constant of dielectric region is set to 4.5so similar to silicon. Our study of the spin-resolved electronic 
structure of 4ZGNR in electric field is based on ab initio local spin density approximation'''^ (ab initio LSD A) and 
norm-conserving pseudo-potentials. Transport properties of the logic gate are studied with non-equilibrium 
green's function (NEGF) technique using Double Zeta Polarized (DZP) basis and the Poisson equation is solved 
self-consistently with integrating the gate-bias into HamUtonian as part of pseudo-potential function™. Transport 
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Figure 1 | Spintronic logic gate of ZGNR. Diagram of the spintronic logic gate and electronic structure of a molecule analogous to the central device of 
the logic gate, (a), The 4ZGNR logic gate consists of four metal GflteS(i,2,3,4) (yellow) and two dielectric regions (red) beside the edges of the ribbon that 
has three parts: left electrode, device and right electrode. Because electrodes are semi-infinite long, only one primitive cell (green) is shown for each 
electrode in this figure. The width for metal gates and dielectric regions is 4 A. Distance between GatCi and Gate2 is 12 A. a is the lattice constant of ZGNR, 
and 2.46 A is used in this letter, (b). The spin density is plotted for the molecule of central device saturated by hydrogen atoms. The isovalue of the spin 
density is set to 0.0358 ;ib/A-' (^b is Bohr magneton), a and P spins are shown in blue and red, respectively, (c). In electric field of 0.0 V/A, 0.5 V/A and 
1.0 V/A, energy levels of the molecule is plotted. Highest Occupied Molecular Orbital (HOMO) and Lowest Unoccupied Molecular Orbital (LUMO) are 
denoted by red and green dash lines. The energy gap between HOMO and LUMO is marked by black arrows. The energy gap of 0.5 V/A is the smallest 
(60 meV) while the energy gaps of 0.0 V/A and 1.0 V/A are around 0.37 eV. 



calculations are carried out in ATK package with k-point sampling of 
1 X 1 X 100 and the cutoff energy is set to 75 Hartree. For molecule 
device of 4ZGNR saturated by hydrogen atoms (Fig. lb), the ground 
state is calculated by DMol' package including all electrons based on 
DNP (Double Numerical plus polarization) basis and LSDA func- 
tional. Based on the first principle calculations of DMoP, the energy 
difference between the ferromagnetic ground state (-1976.518668 
Hartree) and the anti-ferromagnetic ground state (-1976.519288 
Hartree) is 0.62 X 10"'' Hartree (16.9 meV), and the anti-ferromag- 
netic state is confirmed as the ground state. As shown in Fig. lb the 
molecule device is ferromagnetic ordered at each edge and anti- 
ferromagnetic coupled for the ground state. When there is a homo- 
geneous in-plane external electric field (Fig. Ic), the energy gap 
between highest occupied molecular orbital (HOMO) and lowest 
unoccupied molecular orbital (LUMO) is altered and the critical field 
is 0.5 V/A for which the HUMO/LUMO gap is around 60 meV and 
the device is near-conducting. Our result of 4ZGNR molecular device 
is in consistent with previous studies'* ''. When four metal gates are 
buUt up, magnitudes of external electric field at both edges of 4ZGNR 
can be controlled and density distribution of opposite spins can be 
modulated (Fig. 2). Non-equivalent localization resulted from 
applied in-plane bias raises distinct blockade for opposite spins 
and leads to asymmetric transmission spectra. In our prototype, 
input and output signal of the logic gate are defined as bias on metal 
gate and spin-polarized transmission of electrons from left electrode 
to right electrode, respectively. The distance between Gatei and 
Gate2 is 12 A. If a 12 V bias is applied on one metal gate, then the 
electric field nearby this gate is around 0.5 V/A which tunes the 



device into the near- conducting state, whereas the field becomes a 
little larger than 0 V/A when zero bias is applied because the metal 
gates interact with the 4ZGNR and an effective electric field is 
induced due to redistribution of electrons. For input, 12 V and 0 V 
are used to represent 1 and 0 of Boolean algebra. It is expected that 
specific voltages will be used for different devices. 

Gate configurations of input bias are denoted as {yi,V2)^3.^4} (^i 
is the voltage of GatCi in unit of Volt, i= 1,2,3,4). There is a trans- 
mission selection rule and this rule is related to the wavefunction 
symmetry. There are 16 (2'') gate configurations in total. Only seven 
of them are irreducible due to both the mirror symmetries of atomic 
structure along axis and the time-reversal symmetry of spin. For 
symmetric configurations, outputs are the same. Seven irreducible 
configurations are: {0,0,0,0}, {0,0,0,12}, {0,0,12,12}, {0,12,0,12}, 
{0,12,12,0}, {0,12,12,12}, {12,12,12,12} (Fig. 2 a-g). Conducting states 
locate at edges, so only channels at edge are effective for transmission 
around £p. The metal gates next to the edges change the boundary 
conditions of the Poisson equation for electrostatic potential of 
Hamiltonian, and the electrostatic potential changes the effective 
electric field in both the electrodes and the central region. It is the 
inhomogenous effective electric field that changes the distribution of 
spin density in the transport system. For configuration of {0,0,0,0}, 
both forward transmission and backscattering are weak because both 
top and bottom edges are semiconducting, and transmissions of 
opposite spins are both nearly zero. The two transmission spectra 
of opposite spins are non- degenerate because there is a nonzero 
effective field between metal gates and 4ZGNR even under zero bias. 
For configuration of {0,0,0,12}, the top edge is semiconducting and 
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Figure 2 | Transport properties of the spintronic logic gate. Columns from left to right refer to spin-polarized transmission spectra in unit of e^/h, 
external potential of metal gates, transmission pathway and spin density. In the first column, transmission of a and P spins are denoted as grey up-triangle 
and dark down-triangle, respectively. In the second column, the bias of metal gates are displayed and color map shows the corresponding electrostatic 
potential with respect to the color legend from 0 V to 12 V. In the third column, blue arrows denote forward transmission from left electrode to right 
electrode, whereas red arrows denote backscattering. In the fourth column, density of a and p spins are displayed in blue and red, respectively. Rows from 
a to e refer to the seven irreducible combinations of bias on metal gates. 
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Table 1 | Input and output of the spintronic logic gate. Bias on metal Gafei (2,3,41 shown in Fig. 1 a isdenoted as Vi|2,3,4) in unitofVolt. Tajp) 
is the transmission of a(P) spin electrons propagating from left electrode to right electrode, fri+p is the total transmission, /.is the logic value 
of 7" relative to a reference fref- 1 and 0 of Boolean algebra are corresponding to T> T^ef and T< fref, respectively 
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Table 2 | Boolean algebra of the spintronic logic gate: 'AND', 'OR', 'NOT', ip (7ref=0- 1 ) in Table 1 is used as the output signal of the 
device. For inputs (denoted as L;„), 1 and 0 of Boolean algebra are corresponding to 1 2 V and 0 V on metal gates. For outputs (denoted 
as iout), 1 and 0 of Boolean algebra are corresponding to transmission of P spin greater than and less than 0. 1 . Combinations of bias on 
metal gates form all base logic functions in the same device 

.AND. .OR. .NOT. 

{Vini =V, =V2,Vi„2=V3=V4} {Vi„, =V, =V3,Vin2=V2=V4} {V, =V3= 1 2,V2=0,Vi„=V4} 
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transmission along the bottom edge is blockaded due to varying 
chemical potential from left to right which results in a localized spin 
density for p. For configuration of {0,0,12,12}, both forward trans- 
mission and backscattering are blockaded due to the mutative chem- 
ical potential from left to right. For configuration of {0,12,0,12}, the 
top edge is semiconducting whereas the bottom edge is close to 
metallic, so there is one effective channel at bottom edge and the 
transmission of a spin is around 0.3 whereas the transmission of P 
spin is greater than 0.5 because it is close to the resonant peak of 
transmission at energy of — 0.03 eV. For configuration of {0,12,12,0}, 
both edges are backscattered. For configuration of {0,12,12,12}, the 
top edge is backscattered and the bottom edge is nearly metallic. 
However, the backscattering of P spin is much heavier than that of 
a spin due to its very localized spin density, so the transmission of P 
spin is small. For configuration of {12,12,12,12}, channels at both 
edges are open because the localization of spin density is weak, and as 
a result the transmission is big for both a and p. 

For output, a transmission is set as reference T^f, so larger and 
smaller transmission are corresponding to 1 and 0 of Boolean alge- 
bra, respectively. Both spin-polarized and total transmission can 
function in specific logic operations which depend on the reference 
transmission. For example, the logic operation of T^j+p is different 
between T^^{=Q.l and ri.ef=0.6 (Table 1). There is also a remarkable 
difference in output for opposite spins. In {0,12,12,12}, logic value of 
output is 1 for QC spin whereas it is 0 for P spin. Thus, all base logic 
functions are built up in the same device (Table 2) and opposite spins 
can perform separate logic operations simultaneously. Except for 
{0,0,0,12} and {0,12,12,0}, the other five irreducible configurations 
are used for logic operations in the prototype and Tp is preferable for 
output of the logic gate because it has the highest on/off ratio (~7) in 
'NOT' operation. 

In summary, we present a prototype of spintronic logic gates 
consisting of 4ZGNR and four metal gates in this letter, and the 
prototype can function in all base logic operations when specific 
combinations of bias are applied on metal gates. The mechanism lies 



in that, at each edge of 4ZGNR, the ferromagnetic ordered edge state 
around £p can be controlled and modulated by a localized external 
electric field. Combinations of bias make the 4ZGNR running spin- 
tronic calculations as a logic gate. For wider ZGNR, the critical 
external electric field is smaller, thus quantum logic gates of wider 
ZGNR will use voltages lower than 12 V as their inputs and decrease 
the energy consuming. For our spintronic logic gate, the inputs are 
represented by voltages and the outputs are represented by transmis- 
sion coefficients. The transmission coefficients are corresponding to 
the possibility of successful electron transmission through the device 
from the left electrode to the right electrode and thus the coefficients 
represent the transmission current of quantum transport. To make a 
real integrated circuit, there are two methods to convert the output 
signals to input signals. The first method is to use a classical current- 
voltage converting device such as a current transducer to transduce 
the current signal to voltage signal as inputs of other devices. The 
second method is to build up a spintronic current- voltage transducer 
which would distinguish signals of spin-current directly. Based on 
the prototype in this letter, it is possible to design an integrated 
quantum device making more advanced calculations than Boolean 
algebra, like quantum CPUs. 
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